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Based on an equations- of -motion approach for time-dependent pair corre- 
lations in strongly interacting Fermi liquids, we have developed a theory for 
describing the excitation spectrum of these systems. Compared to the known 
"correlated^' random-phase approximation ( CRPA ), our approach has the 
following properties: i) The CRPA is reproduced when pair fluctuations are 
neglected, ii) The first two energy-weighted sumrules are fulfilled implying 
a correct static structure. Hi) No ad-hoc assumptions for the effective mass 
are needed to reproduce the experimental dispersion of the roton in ^He. iv) 
The density response function displays a novel form, arising from vertex cor- 
rections in the proper polarisation. Our theory is presented here with special 
emphasis on this latter point. We have also extended the approach to the 
single particle self-energy and included pair fluctuations in the same way. 
The theory provides a diagrammatic superset of the familiar GW approxi- 
mation. It aims at a consistent calculation of single particle excitations with 
an accuracy that has previously only been achieved for impurities in Bose 
liquids. 



PACS numbers: 67.55.-s, 67.55.Jd, 671.10Ca 



1. MOTIVATION 

Whereas the ground state properties of strongly interacting fermionic 
fluids are nowadays wefl understood, the features of higher-lying excitations 
such as the roton in He and the plasmon in electron liquids still require 
further clarification. Without a proper consideration of pair excitations, 
both the dispersion and the damping of these collective modes are described 
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incorrectly. In the bosonic case, in particular for the roton in ^He, the in- 
troduction of "backflow" correlationJ^ (equivalent to a special case of pair 
excitations) has led to a theor}!^ yielding remarkable agreement with ex- 
periment. The correction of the roton position is due to pair fluctuations 
which have wavelengths comparable to the interparticle distance. 

In Fermi liquids the density-density response function is conventionally 
taken to be of the CRPA form 

X^RPA(<Z,a;) = Xo(g,^)/(l-l>p-h(9)xo(g,c^)) , (1) 

where xo is the Lindhard function and Vp.\i a suitable effective interaction 
(also termed „pseudopotentiai^" and ,, local field^ corrected potential"). Fix- 
ing Vp_-^ through the lj^ and lj^ sum rules leads to a collective mode that 
is energetically much higher than found in experiments. One can cure this 
problem in an ad-hoc manner by introducing an average effective masJ^m* in 
Xo- This is unsatisfactory from the point of view of developing a manifestly 
microscopic description of the low-temperature properties of '^He, and it also 
introduces inconsistencies since it violates the above sum rules. In addition, 
one would expect that the physical mechanisms that lead to the lowering of 
the phonon-roton spectrum in ^He from the Feynman dispersion law are also 
at work in ^He. As in ^He the phonon-roton spectrum has been understood 
quantitatively without the violation of the two above sumrules, our goal is 
here to develop a theory for '^He that has the same level of accuracy and 
consistenty as the '^He theory. 

The need for such a theory is made even clearer in quasi-2-dimensional 
^He, where the collective mode is found inside the particle-hole continuunP. 
An effective m* in an RPA formula may lower (both!) the continuum 
and the collective mode, but only dynamic correlations could reproduce the 
broadening of the phonon. 

2. FORMALISM 

We describe the excited state wave function of strongly interacting 
fermions by a generalized time-dependent Hartree-Fock form 

= ^e-^"*/^Fe^W|<I.o), (2) 

where Af denotes the normalization integral, |<I>o) the Slater determinant 
and Hoo the energy of the correlatecP ground state F|<I>o). The excitation 
operator U includes both single-pair ("c") and two-pair ("d") amplitudes: 

^(0 = J2'^phit)"'p"'h + h dpp'hh'it)alal,af^,a^ , (3) 

ph phh'p' 
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The amplitudes Cph{t) and dpp/hh'it) are determined by minimizing the action 
integral corresponding to the time-dependent Schrodinger equation. Treat- 
ing the external field h^^^ and consequently Cph and dpp'hh' as first order 
perturbations, the resulting equations of motion (EOMs) read 



^^Ph{t) = h;f{t) + {ph'\V^_i,\hp'lcp,h,+ (4) 

{pp'\Vp-h \hh% c;,^,, + {ph'\w\p"p'i dp>'p>,hh' - {h"h'\w\hp'i dpp',h"h' , 

with Cph = ep — Ch denoting the particle-hole excitation energies and W the 
effective interaction in the channels other than particle-hole particle- 
particle (Fp_p), and hole-hole {V^^_^^)■, 

d 



^ — Gph — ep'h' 



dpp' ^hh' 

{pp'\W\hp") Cpf'h' + {pp'\W\p"h') Cpf'h - {ph"\W\hh') Cp>h" - {h"p'\W\hh') Cph" 



+ 



(/iVlV'p-h \p"hX dppn^hh" - {h"p'\V^_^ dppn^h"h' + [IZl 



+ (W'IK-P \P"P"') dp"p"'M' + {h"h"'\V^_^\hh') dpp,^h"h"' ■ (5) 

To simplify the calculation, we retain only those terms that survive in the 
limit of a boson theory: This implies the neglect of exchange effects and also 
of all contributions of the "ladder" type (which are invoked through the Jas- 
trow correlation operator F and are exactly accounted for in the Bose case). 
In the language of commonly used Feynman diagrams this amounts to keep- 
ing the proper polarization parts shown in Fig. ^ Clearly, the graphs no. 3 
and 5—7 invoke vertex corrections, which are beyond the CRPA formalism. 
We also assume a local pair excitation operator, i.e. 

dpp'hh' <4><^l'ah''^h ' Y^d{ci,al) [/5q/5q/ - Pq+q/] , (6) 

qq' 

where /)q is the density operator, and further simplify the c-d coupling by 
replacing the amplitudes Cph{t) by their Fermi-sea averages. With these 
approximations the EOMs can be solved for Cph and dph^p'h' which, in turn, 
yield the linear response wavefunction and thus the induced density 6p{r, t). 
Introducing 

(where = 9{kp — k) = 1 — Uk), the final result for the density-density 
response function can be expressed as follows: 

x''-'i,,u^) = —^^^^^^-- (8) 
1 - K{q,uj)Vp_Y,{q) 



H. M. Bohm, E. Krotscheck, and M. Panholzer 




Fig. 1. Typical proper polarization Feynman diagrams taken into account 
in our approach. (The single and double wavy lines denote the bare and 
screened interactions, respectively. The lines with arrows denote free single 
particle propagators). 



K{q,co) = \- + ^ , (9) 

^ ^ 2^ T?i^ + e(q - q') + e(9 ) 

Here W3 is a three- "phonon" vertex, and e{q) = h^q^ /{2mS{q)) is the Feyn- 
man spectrum. 

We conclude this section by noting that it is easily proved that x{Qt^) 
obeys the lo^ and lo^ sum rules, independent of the specifics of W{q,uj). 
This means that Vp_\^{q) is uniquely defined by the static structure function 
through the first two energy weighted sumrules. 



3. RESULTS 

In Fig. [51 we show the results for the dynamic structure factor S{k, uj) 
of '^He. Input to all our calculations was the static S{k) obtained in Ref. IHl 
The "correlated RPA" (= omitting pair excitations) is shown on the left; 
as mentioned in the introduction the collective mode is higher than found 
experimentally. The right picture gives the CBF result (= including pair 
fluctuations). Obviously, the dispersion is significantly improved and a weak 
pair-excitation background can be seen outside the patricle-hole continuum. 
As the calculated zero-sound mode is practically identical with the experi- 
mental one, a modification of the effective mass is unnecessary to obtain this 
result. 
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Fig. 2. Saturation density ^He dynamic structure factor. Dark line: theo- 
retical zero-sound mode. Filled [open] dots: experimental phonon [magnon] 
modes (Ref. 1^. Squares: RPA. 

It appears interesting to investigate the effective mass of quasi-particle 
excitations in the system within this same formalism and to compare the 
result with the phenomenological m* used in Ref. \^ This can be done 
by studying the time dependent wave function corresponding to a definite 
single particle-hole amplitude, i.e. by omitting the sum over ph (only) in the 
Cpft— contribution in Eq. 121 Then one can proceed in exactly the same way 
as above, leading to a consistent way of studying collective excitations and 
quasi-particle properties. Work in this direction is under investigation. 

In summary, the introduction of dynamic multi-pair correlations, so 
successful in ^He, also proves necessary for obtaining the correct phonon 
dispersion in ^He. A promising extension is the inclusion of spin fluctuations, 
which holds the potential of understanding the large effective mass in ^He. 

As a matter of course, the theory can also be applied to other strongly 
interacting Fermi fluids and to electrons. 
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